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Abstract 

The action of the long-range residual force on the the expectation value of observables in the 
nuclear ground-states is evaluated by finding optimal values for the coefficients of the canonical 
transformation which connects the phonon vacuum state with the (quasi-)particle ground-state. 
After estimating the improvements over the predictions of the independent particle approximation 
we compare the ground-state wave functions obtained using the presented approach with those 
obtained using the conventional random phase approximation (RPA) and its extended version. 
The problem with overbinding of the nuclear ground state calculated using the RPA is shown to be 
removed if one sticks to the prescriptions of the present approach. The reason being that the latter 
conforms to the original variational formulation. Calculations are performed within the two-level 
Lipkin model in which we present results for the binding energies. 
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I. INTRODUCTION 



Approximating the correlated nuclear ground state has been receiving considerable at- 
tention since the early days of nuclear structure physics pQ and still represents a formidable 
challenge. This is an arduous task within the "beyond the mean-field" theories because of 
the action of the residual interaction which brings particle-hole admixtures into the ground 
states. In the present paper we focus on the effects of long-range part of this interaction. 
Unambiguously attributable to the latter is the quadrupole correlation energy, which, as 
shown by the findings in Ref. j2], is considerable and varies between lOOkeV and 5.5 MeV 
in different nuclei. The short-range residual forces on the other hand compete with the 
long-range ones in dominating the ground state shapes formation [3J. As a result of this 
competition in the beginning and the end of major shells the nucleons are paired giving 
rise to spherical shape while in the middle of the shell the nucleons are paired-off and they 
align to the field generating forces thus contributing to deformation. We conjecture that the 
present study can serve as a foundation to investigate the mechanism of the transition be- 
tween these two regimes and in particular on the pairs decoupling process. The mechanism 
that we surmise was concocted in Ref. |3j and essentially implies that the long-range force 
breaks nucleon pairs which may further recouple due to the pairing force. 

We approximate the nuclear ground state wave function with the phonon operators [5] 
vacuum state. A general form of the phonon vacuum was proposed by Sorensen [6] and later 
Goswami and Pal [7] estimated explicitly the correlation coefficients of the 2p-2h admixtures 
into the BCS wave function [S] relating them to the forward and backward phonon ampli- 
tudes. The expression they derived turned out to be also valid for higher order correlations 
[9] in the the random phase approximation (RPA) pU]. Being a small amplitude limit of the 
time- dependent Hartree-Fock approximation, however, RPA is known to be able to account 
only for small correlation effects. Since in open-shell nuclei the backward phonon amplitudes 
are by no means small, RPA is becoming questionable in describing the low-energy states 
of such nuclei. This problem was addressed by Hara jU [T2] who proposed an improvement 
over the RPA based on the Pauli blocking principle which plays a progressively important 
role with the increasing number of the valence nucleons. This extension has proven to be 
in better accord with the experimental data as demonstrated in a number of papers as for 
example in Refs.[13l [H]. Although its superiority over the standard RPA is undeniable, 
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the variational character of the theory is violated as the ground state is found to be over- 
bound. The strong argument that the variational property of a theory insures a converging 
succession of approximations to the exact solution fostered the formulation of a elaborate 
formalism, called self-consistent RPA [15) . which as in the conventional particle- hole theory, 
allowed to take into account the nucleon correlations without explicitly constructing the 
ground state wave function. 

In the present paper we keep using the explicit form of the fermionic many body vacuum 
[0] but depart from varying the excited state wave function. On the contrary we use the 
correlation coefficients as parameters which we fix by optimizing the ground state trial wave 
function using a variational procedure. This approach benefits from the findings in Ref . [16] 
where it was shown that this class of wave functions is a vacuum for a generalized phonon 
operator, adding to the standard one specific two-body operators correcting for the Pauli 
principle. By way of example using the two-level Lipkin-Meshkov-Glick(LMG) model [17] 
they showed that the additional terms improve the convergence substantially. In this way 
the phonon vacuum state absorbs additional correlations effects than of ones obtained using 
the equations-of-motion method for the standard phonon operator. 

The paper proceeds as follows. In Section [II] we outline the problem and summarize the 
main obstacles towards the exact solution. Basic equations of several approximate methods 
including the RPA, ERPA and the explicit variation of the phonon vacuum state along with 



the exact solution within the LMG model are derived in Section III A comparison between 
them is performed on the basis of the ground state total energy in Section IV Summary 
and outlook is given in Section |Vj 

II. FORMULATION 

Formally a wave function which contains quasiparticle admixtures into the independent 
quasiparticle wave function and is a vacuum for the phonon operators [5] 



= \ E {^m'Wi - ( n A 'vi;-n././'A _ ^ ( i) 

a' 

can be expressed as [6] 

I) = N J) (2) 
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with 



S = -\J2 SxUiJ^j'^-f^A^j^X^U^X - n). (3) 

12;A/t 

The coefficients S\(jij[ \j2f2) are referred to as correlation coefficients and denote the 
amplitudes for the presence of zero, four, eight, ... quasiparticles in the ground state due 
to the virtual vibrations. These present a primary source of structure information for the 
ground states and make up a major part of our present research. The correlated and un- 
correlated ground states are denoted as |) and ) correspondingly. The quasiparticle-phonon 
nomenclature follows the one given in Ref. [5]. 

Using the ERPA, the correlation coefficients are found to satisfy the equations [U E] 

= E(* - PhiO^OdiM)^ (4) 
where pj is the quasiparticle occupation number on the level j. 

In Eq.@ iV is a normalization factor which in physical terms is the overlap between the 
independent-particle and the correlated wave functions. It is found to be 

N ° = ( e (st+s)) ( 5 ) 
In RPA, suggesting small correlations so that higher order terms contribute relatively 
little, this constant is approximated as : 

TV 2 ^ - (fil 

e lE J1 ^ ;A (2A + i)s^(MU 2 ii)- 1 } 

An explicit solution to the system Q was obtained by Hara jl]. 

Changing the frame of mind we shall try to obtain the correlation coefficients explicitly 
varying the wave function |) in the functional 

m\) = o (7) 

, with <Sx (j 1 J i I ) being variational parameters, i.e. we shall try to solve the equation 
5(N*(e st He s )) = (e s ' He s ) (6N*) + iV 2 (<5(e 5t He s )^j = 0, (8) 



4 



with respect to 5*a (jiJi I J2J2) - If one restricts the configuration space to 2p-2h admixtures 
only, the quantities that need to be evaluated are presented in the expression below 

st s (H) + 2(HS) + (ggg) 

Evaluating these quantities proves prohibitively laborious in the realistic case even within 
the restricted space and some approximate expressions for these will be given elsewhere. In 
order to assess the utility of our approach we have recourse to the widely used LMG model. 



III. SOLUTION WITHIN THE LIPKIN-MESHKOV-GLICK MODEL 

In order to access the utility of different approaches and to prove the usefulness of the 
proposed scheme we limit the configuration space and simplify the inter-nucleon interaction 
to monopole-monopole one as suggested by Lipkin, Meshkov and Glick [TTJIIH]. This setting 
permits comparisons between the rates of convergence of different approximation methods, 
including the hereby described, to the exact solution. 

The interaction of N particles on 2 quantum levels is presented by the following Hamil- 
tonian 

H = H + V; H = eJ ; V = ^(J + + J_) 2 (10) 

where 

J+ = a\ a a>-ia, J- = o}_ la ai a , Jq = - '^ j (o\ a o,ia — d-iofl-i<x) (11) 

a a a 

are analogous to the raising, lowering and angular momentum' ^-component of the quasi- 
spin algebra respectively, a) represents the particle creation operator, ±1 in the subscript 
denote the upper or lower level, e is the energy gap between the two levels and G is the 
interaction strength. 

We shall also make use of the operators 

s t = 4a a -ia; s~ = al la a la ; s° a = ^(a\ a a la - al la a_ la ). (12) 
The exact solution 
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|¥>=5>|* 4 ) (13) 

i 

is a superposition of states with 0, 1, 2 . . . , iV particles on the upper level (|$o) =))■ The 
weights Cj are readily obtained by solving the eigenvalue problem 

J2(®n\H\$ n >)c n , = Ec n . (14) 
n' 

The non-zero elements of the matrix on the left-hand side of the above equation evaluate 



to 



(K\H\* n ) = ( -j + n}e + G(~ + n + (n + l){N-n) ) (15) 



G 



($;|#|$ n+2 > = -y/n(n - 1)(N + 2 - n)(N + 1 - n). (16) 
For the ground state total energy we then obtain 



E = ^($Jff* < ) + 2^c^ J (tJff* w )- (17) 

i i 

We shall further present the solutions for exited states containing only one particle on 
the upper level and one hole on the lower one, i.e. 

\lplh) m = (X>- m)s - - ^ m) ^)) ) ( 18 ) 

Employing the RPA, i.e. 



{K,s+\)=8 aa , (19) 

one obtains the well-known equation for the excitation energies and the forward and 
backward amplitudes: 



8 aP {e - G) + G -Gil-Sap) \ hp \ I \ ip 



id 



-G(l-M S afj (e-G) + G \<p \0-I \p 



(20) 



which together with the normalization of the wave functions (18) yields 
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(0) 



1 1 + x/2 + coo/e 



X/2 



(21) 



where 



2G(iV - 1) /- 

X = > c;o = £V 1 + X- 

The ground state total energy evaluates to 

E = E HF - Nu {if ) 2 , 

where 



(22) 



(23) 



E 



HF 



(G-e 



N 



(24) 



The condition (19) disregards some aspects of the nature of the excited states (18), in 
particular the fact that the number of particle-hole state in the ground state may be non- 
negligible if sufficiently strong interaction is applied. In a broader context, than the hereby 
considered, Hara [I] suggested to include explicitly the number of quasiparticles on each 
level, which turned out to have dramatic effect on the collective properties of the low-lying 
states in even-even nuclei [321 EH]. Adapting it to the LMG model we can write 



(\n +1 \) = Np, (\n^\) = N(l-p). 



Eq. ( 19 ) then transforms to 



(25) 



|s-,s+i> = <U'(l-2/>). 



Analogous to the equation (20) in the current context is the following one 



A R 

TJ* A* 











H 













-TJ. 



OL\OL2 




where 



(26) 



(27) 



A aia2 = Gil - 2pf - 5 12 (1 - 2p)(G - e), (28) 
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and 



The solution is 



B aia2 = G(l - 2p)(5 aia2 + 2p- 1), (29) 
U aia2 = 6 aiaa (l - 2p). (30) 



^(o) = 1 l + x/2 + o;o/e 

^= - 1 -= X/2 (32) 
v /iV(l-2p) v /( 1 + f ) + x + 



where 



of 

Wo 2 = e 2 (l + X ), x = _[(l_2p)iV-l]. (33) 
The system of equations closure is insured by adding the following one: 

(v°f ,,,, 

p= rT%^- (34) 



Finally, the featured method that we examine (conf. Sec. [IT]) translates in the language 
of the LMG model in the following way. The wave function ^ assumes the form 

|) = N e^ a i a i Saia2S ™ lS ^}. (35) 

If we truncate the exponent expansion to second order, i.e. allow for 2p-2h admixtures 
only into the ground-state wave function, we obtain 

^ = TTTv rhr* — ToT" • ( 36 ) 

The variational problem then is rewritten as 

5 (n 2 /(l + i J2 S aia2 s ar s ai )H(l + l -Y, S ai aAA 2 ))) = 0- (37) 

\ \ 0:10:2 aict2 / / 

It can be rigorously proved, that in this simple model the correlation coefficients are 



all equal, i.e. 5y = S. Then performing the variation (37) leads to the following simple 
quadratic equation for S: 
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1+ (2^ + 2N S -^(N 2 - N)S 2 = 0. (38) 
The ground-state energy in this case evaluates to 



E = N 2 



E HF + N(N - 1)GS + J(e - G)(-N + 4)N(N - 1)S 2 + GN(N - 1) 2 S 2 



• (39) 



Performing the variation ([7| analytically in an extended configuration space represents 
a formidable challenge even in this simple model. With the aid of computers, however, we 
managed to evaluate the quantities of interest in a 4p-4h space for systems with up to 8 
particles as will be shown in the last section. 

IV. NUMERICAL COMPARISON 

The three approximations presented in the previous section are compared with the exact 
solution based on the binding energies (FigjIJ. From this figure one clearly designates the 
critical RPA strength 

Gcrit = ~2(N-iy (40) 

which delineates the regions where a real solution exists and the one in which only a 
complex solution can be obtained. It is worthwhile to notice that the critical point stands 
just at the transition between the two nearly linear sections of the exact solution which are 
more distinguished in systems with a larger number of particles. The ERPA solution is also 
plotted in the same region as the RPA, i.e. (G cr it, 0] and there ERPA performs much better 
than RPA as the former does not exhibit the overbinding near the critical point of the latter. 
Formally one can derive a critical strength for the ERPA as well: 

Gcrit = ~2[(l-2p)N-iy (41) 
The simplicity of this formula, however, can be deceptive because the particle occupation 

p depends on the interaction strength, i.e. p = p(G), which implies the self-consistent nature 

of the problem. We left this problem open for a future study. 

The explicit variation of the phonon vacuum with 2p-2h admixtures on the other side 

yields solution for any G. Increasing the interaction beyond the RPA critical point causes 
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FIG. 1. Binding energies (MeV) of two-level LMG model systems with N = 8, 20, 50, 100 particles 
as function of the interaction strength G(MeV). The energy gap between the two levels is set to 1 
MeV. 

progressive divergence of this solution from the exact one. This divergence exacerbates 
incrementing the number of particles in the system which elucidates the fact that in the 
strong interaction regime multi-particle-hole admixtures start to play a dominant role. Our 
endevours to take into account 4p-4h configurations have been impeded by the calculation 
of the following average values 

(s 2 SlS3S4J-J+4444)> (42) 

which we managed to handle with the aid of a computer for systems of up to 8 particles. 
The result is presented in Fig [2] which exhibits the fact that the 4p-4h solution is closer to 
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- Exact 

- - Phonon vacuum 2p2h 

- - Phonon vacuum 4p4h 

-0.16 -0.14 -0.12 -0.10 -0.08 -0.06 -0.04 -0.02 0.00 

Interaction strength 

FIG. 2. Same as Fig. [I] but for N = 8 particles only with 4p4h configurations included. 

the exact one than the 2p-2h solution which vindicates the prescription of the variational 
principle for monotonic convergence. 

V. CONCLUSION AND OUTLOOK 

In this work we initiated the development of a variational approach for approximating 
the ground state of nuclei tailored to take into account the action of the long-range residual 
forces. Essentially it is an attempt to provide a controlled succession of approximations for 
estimating the contributions from different multi-particle-hole admixtures. Applying our 
idea to the uncomplicated LMG model we showed its superiority over the RPA based on 
comparison with the exact solution for the binding energy. One drawback of the present 
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development is that the complexity of calculating analytically higher order terms, which 



as shown in Sec. IV play increasing important role in systems with larger number (> 10) 
of particles, escalates rapidly. For that purpose specific summation methods might prove 
very helpful. Being a many-body wave function the present development is expected to 
excel in describing phenomena involving many-particle correlations as, for example, cluster 
configurations. The influence of the long-range residual forces on the mean-field and the 
pairing correlations is another perspective which is currently progressing. 
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